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Abstract 

The differences between subcritical liquid drop and 
supercritical fluid drop behavior are discussed. Under 
subcritical, evaporative high emission rate conditions, 
a film layer is present in the inner part of the drop 
surface which contributes to the unique determination 
of the boundary conditions; it is this film layer which 
contributes to the solution’s convective-diflFusive char- 
acter. In contrast, under supercritical conditions the 
boundary conditions contain a degree of arbitrariness 
due to the absence of a surface, and the solution has 
then a purely diffusive character. Results from simu- 
lations of a free fluid drop under no-gravity conditions 
are compared to microgravity experimental data from 
suspended, large drop experiments at high, low and in- 
termediary temperatures and in a range of pressures en- 
compassing the sub- and supercritical regime. Despite 
the difference between the conditions of the simulations 
and the experiments (suspension vs. free floating), the 
time rate of variation of the drop diameter square is re- 
markably well predicted in the linear curve regime. The 
drop diameter is determined in the simulations from the 
location of the maximum density gradient, and agrees 
well with the data. It is also shown that the classical 
calculation of the Lewis number gives qualitatively er- 
roneous results at supercritical conditions, but that an 
effective Lewis number previously defined gives quali- 
tatively correct estimates of the length scales for heat 
and mass transfer at all pressures. 

Introduction 

Supercritical fluids are a common occurrence in propul- 
sion systems, yet the fundamentals of their behavior is 
poorly understood. This is due to the unusual proper- 
ties of supercritical fluids which combine the character- 
istics of both liquids and gases, making their behavior 
unintuitive. Furthermore, data obtained from (mostly 
optical) measurements can be misleading because re- 
gions of high density thus observed are frequently iden- 
tified with liquids: a common misconception is that if 
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in an experiment one can optically identify ‘drops’ and 
‘ligaments’, the observed fluid must be in a liquid state. 
For example, Mayer et al. [Ij and Krtlile and Mayer [2] 
concluded that the LOx was in a liquid state. This 
inference is generally incorrect because in fact optical 
measurements detect any large chfinge (i.e, gradients) 
in density. Thus, the density ratio may be well below 
O(IO^) that characterizes its liquid/gas value, but the 
measurement will still identify a change in the index of 
refraction providing that the change is sudden (steep 
gradients). As shown by simulations of supercritical 
fluids [3], the density gradients may remain large dur- 
ing the supercritical binary fluids mixing, thus mak- 
ing them optically identifiable. Therefore, there is no 
inconsistency between the optical observation of high 
density regions and the fluids being in a supercritical 
state. A second misconception is that because a fluid 
has a liquid- like density, it is appropriate to model it as 
a liquid. However, such fluids may have liquid- like den- 
sities while their transport properties differ from those 
of a liquid. In fact, it was more recently observed and 
recognized that supercritical fluids behave very differ- 
ently from liquids [4], [5] and [6] in that their breakup no 
longer resembles the well known atmospheric pressure 
atomization transforming a liquid sheet into approod- 
mately spherical drops. 

Considering that the critical pressure of most fuel 
hydrocarbons used in Diesel and gas turbine engines 
is in the range of 1.5 - 3 MPa, and the fact that the 
maximum pressure attained in these engines is about 6 
MPa, it is clear that the fuel in the combustion chamber 
will experience both subcritical and supercritical con- 
ditions. Studies of drop behavior over a wide range of 
pressures were performed in the past (Yang et al. [8], 
Delplanque and Sirignano (9), Haldenwang et al. [10], 
and the review of Givler and Abraham [11]), however 
none of these studies identified the crucial differences 
between the subcritical and supercritical behavior. In 
fact, in two of these studies [8], [10], it was found that 
the subcritical and supercritical behavior is similar as 
the drop diameter decreased according to the classical 
d^-law [12] over a wide range of pressures and drop di- 
ivmeters. Moreover, despite the existence of micrograv- 
ity data for model v 2 ilidation (see below), none of the 
existing models has undertaken a thorough comparison 
between numerical predictions and observations. 



rUo prnsout study is (l(jv(^trd to tlu» exploration of 
ddfernnros in Huid-lxdmvior rharactnristu’s under sub- 
critical and suijercritical conditions, and to the model 
vjilidation in both of those regimes. We choose to exer- 
cise the model for heptane fluid drops in nitrogen be- 
cause of the availability of data for model validation. 

Model equations 

The configuration studied is that of a single spherical 
drop in a medium with specified far field conditions. 
These far field values are identified by the subscript 
‘e’ and the location of the far field boundary, Re{t), is 
calculated in a Lagrangian way to be that of null mass 
flux. 

The conservation equations are based upon Keiz- 
er’s fluctuation theory [13] which has the distinct ad- 
vantage of accounting for non-equilibrium processes. 
This formalism therefore leads to the most general fluid 
equations where the partial molar fluxes, J i, and the 
heat flux, are related to thermodynamic quantities 
as follows: 

N N 

7 . = T = 

J j 

( 1 ) 

where /3 = T is the temperatxu*e, is the 

universal gas constant, and fij are the chemical poten- 
tials. Here Lij are the Pick’s diffusion elements, Lqq is 
the Fourier thermal diffusion element, Liq are the Soret 
diffusion, Lqj are the Dufour diffusion elements, and the 
Onsager relations state that Lij = Lj^ and Liq = Lqi. 
Additionally, conservation of fluxes and mass in the sys- 
tem imply that ^ and Lijrrii = 0 for 

j € [1, iV] and j = q where are the molar masses and 
N is the total niimber of species. 

Using the thermodynamic relationship 

d{0Hj) = Hivjdp - hjd\nT) + ( ^ aDiidXi)IXj (2) 

1 

where is the partial volume, p is the pressure, Xi is 
the molar fraction, hj is the molar enthalpy, and the 
mass diffusion factors are 

aD^J = HXidpJdXj = dX,/dXj^X,d\xi-fJdXj (3) 

one can calculate J » and ^ from 1 and 2. The conser- 
vations equations were derived in detail in Harstad and 
Bellan [3| and results were obtained there for supercrit- 
ical conditions. It was shown there that the transport 
matrix can be calculated as a function of the thermo- 
dynamic variables and of two transport matrices: one 
transport matrix is symmetric and is called the mass 
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mat fix. The iurru aeies of the i)ro[;er dehnition of trans- 
[)ort coetfieitnifs tliat an' coiisistimt with experimentally 
UK'JVsured ones w.us addn\ss<Hl in Harstad and Bellan 
[14]. The emphasis is liere on extending the calcula- 
tions with the same model to the range of subcritical 
conditions. It turns out that one of the major differ- 
ences between subcritical and supercritical regimes is 
in the boundary conditions at the drop boundary; it is 
this condition which is discussed below in detail, the 
two other boundary conditions (at the drop center and 
in the far field) being trivial. 

Boundary conditions 

The detaiiled boundary conditions at r = Rfi have been 
derived in [3|; here we describe only the new aspects 
that enabled the calculations to be extended to sub- 
critical conditions. The jump conditions at the drop 
boundary are: mass balance; relationship between Rd 
and the emission flux Fcms\ heat balance; balance of 
species 1 flux; and the nonequilibrium evaporation law. 
Additional equations at r = are the momentum, 
auid the equation of state for the mixture which is 
used twice (once on each side of the boundary). Thus 
there is a total of eight equations and nine unknowns: 
-^ 16 - ^ 16 - Pfc ' P?- ^d, Tb, Pb and where the sub- 
script b denotes the drop boundary, p is the density, and 
superscripts L and G denote the inside and outside re- 
gions of the drop; variable is obtained by integrating 
the drop continuity equation starting at r = 0. A ninth 
independent relationship exists only under subcritical 
conditions and is related to the existence of a surface, 
as discussed below. 

The indeterminacy of the boundary conditions for 
a fluid drop under supercritical conditions has already 
been discussed by Harstad and Bellam [3] . This is physi- 
cally understandable since there is no true surface, and 
thus there is an arbitrariness as to the choice of the 
boundary to follow. At least three choices are reason- 
able: One may follow the pure fluid boundary as was 
done by Harstad and Bellan [3]. Another possibility is 
to follow the initial boundary separating the two fluids, 
this being the choice in the present calculation. The 
third possibility is to follow the point of maximum den- 
sity gradient; although this is not the present choice, 
the point of maximum density gradient is calculated 
here d posterton to indicate the location of the optically 
identified fluid drop. In contrast, under subcritical con- 
ditions the boundary to follow is the drop surface and 
the problem is fully determined. 

There are other important consequences of the ex- 
istence or lack of a surface at r = /?^. For example, un- 
der strong evaporative conditions a mass fraction “film’ 
layer exists inside the drop [15] and the thickness of 
this layer, 6y Ar“ where Ar” is the distance from 
the surface to the first grid point inside the drop. A 
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of D^ff and that of an e<^uivalent thermal conductivity 
\^ff were calculated under the quasi-steady assump- 
tion in [7] by finding two linear combinations of T and 
the mass fraction V'l for which the transport matrix can 
be approximately diagonalized. In diagonal form, the 
characteristic length scales for diffusional transport of 
these two new variables are apparent, and this allows 
the definition of £)«// and [7]. Previous calcu- 
lations [7] show that Ag// > A and that D^f/ < D. 
These definitions also allow the calculation of an ef- 
fective Lewis number Le^/j = Xgffj{nCpDt,ff) once 
the values of the dependent variables are known. The 
quasi-steady assumption does not remove the general- 
ity of the estimate since the essence of the estimate is 
that of a characteristic length. 

One of the most important consequences of the 
mass fraction him layer existence is the direct rela- 
tionship that exists between Y{Rd — ff) and Y{Rd + e), 
where e max(<5y,<57'); it is this relationship which 
provides the needed additional equation to fully de- 
termine the solution at the drop surface. This re- 
lationship can be formulated by considering the dif- 
ference AY^^ = Yi^{Rd - e) - Yy^{Rd - Ar") where 
Yi^{Rd - Ar“) represents the computational grid cen- 
ter value at the first adjacent position to the film layer 
inside the drop such that > 6y- Similarly one may 
define AT^ = T^{Rti — £)—T^{R^-Ar^). The variable 
= exp{ifj - where is the fugacity, quantifies 
the Yj jump across the drop surface and can be cal- 
culated from the state equation. For example, under 
strict equilibrium evaporation (i.e. Fems = 0) condi- 
tions, = 1. For finite F^ms smd for a binary mixture 


system, its ratio to a reference state Fre/{^ 1 ,^ 2 ) 
defined by 

can be 

F mrns ~ ^pFref 

( 4 ) 

where 



( 5 ) 

Bj = aajrrijUTjVp , j = 1,2 

( 6 ) 

- ^) = (^2 ~ l)/(v 2 - ^i) - 

( 7 ) 

X^ilU +S)= ^,(^2 - l)/(^2 - 

(8) 

~ — (1 ~ Sl)/(^2 ~ ^ l) + 

( 9 ) 

+ -) = ^ 2(1 ~ ^l)/(^2 “ ^l) + ^av^F 

( 10 ) 


where aaj is an accommodation coefficient, utj is the 
mean molecular velocity crossing the plane in one di- 
rection, and n is the number of moles per unit volume 
[3|; consistently, Fre/(l, 1) = 0. A detailed analysis of 
the film layer yields then 

+ ( 11 ) 

^r-)]/[B2Yl‘(RH - Ar-) + - Ar")] 


vvliK'h [irovidi'h flu* atldir ioiuil rc‘lari()u.shi|> tluit allows 
closure' of rlio systt'ui of o<iuaf.ions at the* drop bound- 
ary. Since' imdor supf'rrritical conditions the concept 
of latent he'at, aiul therefore of evaporation, is not aj> 
plicabl(', tli(' above analysis df)es not hold. However, 
the film layer computatiotial approach is still necessary 
if the Peclet numl)er of the grid Pegnd > 0(1) in or- 
der to insure that all scales are resolved. Therefore, 
the formalism of the film layer is retained for computa- 
tional purposes even under supercritical conditions, al- 
though the layer no longer exists physically. Essentially, 
the solution in the supercritical regime has a diffusive 
character, whereas in the subcritical regime it has a 
diffusive-convective character where the convective part 
is introduced by the film layer and the evaporation. 

Results 

The present simulations are performed for an n-heptane 
drop in nitrogen because it is the set of binary sub- 
stances which is best documented experimentally. The 
equations of state have been calculated according to 
the procedure described in Harstad et al. [16], and the 
calculation of properties has been described in Harstad 
and Bellan [3|. The purpose of the present simulations 
is to validate the model; we also present here some re- 
sults emphasizing the importance of the effective Lewis 
number concept. 

The only data that can be used for comparisons 
is that obtained under evaporative rather than burning 
conditions, since in the last case the flame temperature 
that acts as the far field boundary is unknown. Further- 
more, as shown below, it is only microgravity data that 
can be considered valid for these comparisons because 
normal gravity data has unavoidable convective effects 
that are not modeled here. However, Chauveau et al. 

[17] have shown that for the 10^^ - 10“"^ ratios of the 
microgravity to normal gravity obtainable both in par- 
abolic flights and in drop tower experiments, buoyancy 
effects on the evaporation constant can still be iden- 
tified. Additionally, since all high pressure micrograv- 
ity drop evaporation experiments were performed with 
suspended drops, even these data are clearly not totally 
equivalent to our simulation results which are obtained 
for a free floating drop. The recent estimates of Morin 

[18] show that the effects of the suspending fiber in- 
crease with increasing pressure and temperature, and 
with decreasing drop size. 

To our knowledge, microgravity obtained data with 
CtHiq drops evaporating in N 2 were reported only by 
Sato [19] and Nomura et al. [20]. In their experiments 
the 0.7-1 mm drops were suspended from a fiber of at 
least 100;i diameter. The CjHiq drop evaporation ex- 
periments of Chauveau et al. [21] were conducted only 
in normiil gravity, whereas their reported micrograv- 
ity experiments were of burning drops. Therefore, our 
comparison focuses on the data of [19] and [20], while 
i\lso considering for reference (see Table 1) the more re- 
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TIh^ simulation.s won^ porformnd for nominal initial 
t OIK lit ions matdiing f lu' <*x|KTii[iniital data; /?JJ = 0.35 
min nxcopt for tlin comparKson with Sato’s [19] data 
which wiis performed for f?Jj =0.5 mm, and = 300 
K. The far hold conditions are located at = 4 mm 
whore and p*. are specihed consistent with those of 
the experiments and = 0. The fluid drop is initially 
composed of pure heptane (Tc = 540.3 K,pc == 2.76 
MPa), while the surrounding is nitrogen (Tc = 126.2 
K, Pc = 3.39 MPa); in order to avoid an initial unphys- 
ical discontinuity, a minute amount of heptane exists 
initially in the drop surroundings, its distribution van- 
ishing with increasing r. For the same reason, although 
the fluid drop temperature and drof>-surroundings fluid 
composition are assumed initially uniform, a set of com- 
putational initial conditions (i.e. spatial profiles of the 
variables) are calculated for each simulation by satisfy- 
ing the nominal initial conditions at the domain bound- 
aries and the jump conditions at 

In all of the discussions below, ‘subcritical’ and 
‘supercritical' qualifications will be used with respect 
to the heptane critical point, and not with respect to 
the critical point of the mixture which varies according 
to the local composition. 

Determination of thermal diffusion factors from 
high temperature data 

To proceed with the calculation, one must specify val- 
ues for the thermal diffusion factors, a, which can be 
defined either from the Irwing-Kirkwood (IK) or the 
Bearman-Kirkwood (BK) form of the heat flux [14] and 
OiBK.ij = oc[K,ij - ctk,ij where 

= (m^7nj/m)(hi/mi - hj/mj)/{R^T). (12) 

However, values of asK are not well known for most 
substances, except at atmospheric conditions where 
they can be calculated from kinetic theory. Since we 
are here interested in calculations at considerably larger 
pressures, the question arises as how to calculate Qbk^ 
For this purpose, the premise is that if it can be shown 
that o/zc is very small, in fact it can be considered negli- 
gible with respect to aa/c -q/k, auid then aa/c — -cch 
which is only a fimction of thermodynamic quantities 
[14|. Since is calculated from thermodynamics, this 
would provide an approximate value for asK for all 
(p, T) conditions where au^joY is very small having de- 
fined a* = max(^p x,x,) I Of/i |. Heptane/nitrogen plots 
of a* illustrated in [14] show that a* ~ 5 at 0.1 and 0.5 
MPa for T € [300 A', 800 A'] . Therefore, it is expected 
that aix/ct* 1 is indeed correct in the range of (p, T) 
of interest. The purpose of these high temperature data 
comparisons is to verify the consistency of the premise 
that ocfx/ot* ia very small. A similar premise for ocbk 
is also tested. 
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lations portraying Nomura ot al.’s [20| experiments at 
liigh temperature (745 K) in tlu* i>ressure range of 0.1 
- 2 MPa. lu tlie muiKTical simulations, the location of 
tli(' <lro[) [)oun(larv is defined to be that of the maximum 
density gradient, for consistency with optical measure- 
ments. In agre^mumt whth well known theory [12], at 
0. 1 MPa the liquid/giu» interface is found to be precisely 
that of maximum density gradient. With increasing p 
the two locations still coincide for all simulations in the 
range 0.1-5 MPa investigated in this work, but the 
density greidient, although still substantial, decreases 
across the boundary as p increases. 

Ail but two of our simulations were conducted with 
= 0.01; the remaining simulations was conducted 
with aBK = 0.01. Our results with a//c = 0.01 capture 
the 0.1 MPa data very well but display a somewhat 
earlier d^-law behavior; it is unclear whether the non- 
coinciding part of the data and simulations fall within 
the experimental error since this error is not provided 
with the data. However, the results for 0.1 MPa with 
= 0.01 depart considerably from the data. Sim- 
ilarly, the data at 0.5 MPa is compared with results 
from simulations with both asK and a/^ specified as 
0.01. It is clear that the aBK = 0.01 results fall short 
of agreement with the data, and in fact show a typical 
large increase in the evaporation time that was obtained 
with as/<' = 0.01 at other pressures as well. In contrast, 
the ajK = 0.01 results capture the nonlinear portion of 
the curve very well with a small discrepancy in the total 
evaporation time. Simulations and data at 2 MPa agree 
only during the initial time, after which the simulations 
display the expected smooth variation consistent with 
drop heating, whereas the data exhibit two disconti- 
nuities that can be explained only by the presence of 
the suspending fiber. Calculated slopes of the linear 
part of the curves, called the evaporation constant [12], 
AT, are presented for comparison in Table 1 for the 0.1 
and 0.5 MPa results obtained with ajK = 0.01. Despite 
the presence of the suspending fiber in the experiments, 
an excellent agreement exists between simulations and 
data. A similar comparison cannot be performed at 2 
MPa since there is no evidence of linear behavior in the 
data. 

Confirmation of thermal diffusion factors from 
intermediate temperatures data 

Displayed in Fig. 2 are p = 2 MPa comparisons of sim- 
ulation results at 655 K for various values of one 
simulation where asK instead of a/zc is prescribed, and 
Nomura et al.’s [20] data at 656 K. The numerical pre- 
dictions are a very weak function of aiK in the range 
-0.6 - 0.6 and agree remarkably well with the data dur- 
ing the initial heat up period of the drop. Eventually, 
the data shows a faster eviiporation than our simula- 
tions. although the lack of error bars in the data make 
it impossible to evaluate the extent of the disagreement. 
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(during the pxporimont ) on fhe evaporation process. 
H(wever, results with a/y/v - 0.0 1 clearly ovr'restiinate 
both the growtli of the drop diirinj^ the initial heat up 
time and th<' drop <*vaporation time. These results are 
consistent with those of Fig. 1. 

Additional comparisons between numerical predic- 
tions and data is portrayed in Fig. 3 where comparisons 
are made in the range 0.1 - 2 MPa between simulations 
at 655 K with a/K =0.01, and data in the range 648 - 
669 K. The initial heating time is again very well repro- 
duced by the simulations, except that the predictions at 
0. 1 MPa display again an earlier ci^-law behavior. The 
evaporation time is very well reproduced at 0.1 MPa, 
and less well as the pressure increases. Since it is dif- 
ficult to quantify the influence of the suspending fiber 
on this data as the pressure increases, we can qualify 
this comparison as very encouraging. 

Table 1 includes comparisons of K for this interme- 
diary temperature regime, and shows good to excellent 
agreement between data and predictions. 

This study indicates (see also below) that the value 
of aff(/a* is indeed rather small and that asK — 
is correct. The assumption made in all calculations pre- 
sented below is that the value of arn: is the same small 
value determined at high temperatures regardless of the 
(p, T) conditions, and thus that asK — *~a^.This as- 
sumption might not be entirely valid, as in general a//c 
is a function of both p and T. This assumption, the 
fact that the data is from suspended drop experiments 
whereas our calculations are for free drops, and the 
slight but nevertheless present influence of buoyancy 
for the stated microgravity level of 10”^ 10^^ in the 

experiment might explain the 15-20% discrepancies (see 
below and Table 1) betw^een data and results from sim- 
ulations. 

Comparison with data at low temperatures 

The low temperature data of Nomura et al.[20] and Sato 
[19] (Sato’s data was approximated from his figure) is 
shown in Fig. 4 along with our numerical predictions 
at 445 K, 470 K and 495 K using ajx — 0.01. The 
temperature range for Nomura et al.’s [20] data is 466 
- 493 K whereas Sato’s [19] data was obtained at 445 
K; the data in [20] is in the 0.1-5 MPa range, whereas 
that of [19] is at 2 MPa. The comparisons are very 
good at low p and deteriorate as p increases. The pre- 
dictions and data [20] agree remarkably well at 0.1 and 
0.5 MPa, whereas at 1 MPa the evaporation time is 
slightly overpredicted by the simulations. Nevertheless, 
the calculated and measured evaporation constant (Ta- 
ble 1) show very good agreement at all three pressures. 
The 2 MPa numerical results approximate the ex- 
perimental variation [19] fairly well, and the agreement 
ill the value of K (Table 1) is excellent. At p = 5 
MPa, our simulation of a free drop shows an increased 
heating time, whereas the suspended drop in the exper- 


iment slunvs ri < lr< rea.se(l heating time with res(x*(*t U) 
tlie 0. 1 MPa case. Tin' difference betweeii tin' (Experi- 
mental c(;iulif.ions and those of the simuiarions (Explains 
the disagr('(nn«mt in th(' iK'at up tilin', although the rate 
of regression of tho largest gradi('iit location is surpris- 
ingly well predic:t('d. Since at 5 MPa the conditions are 
supercritical, tlu're is no evaporation and the concept 
of evaporation constant is irrelevant, although compar- 
isons between the rates of regression are still meaning- 
ful. 

Determination of the Lewis number 

Figure 5 (a - d) illustrates the classical and the effective 
Lewis number, Le and Le«/ / , for the high temperature 
calculations (745 K) validated above. Figures 5a and 
5b show results from the 0.1 MPa simulations, whereas 
Figs. 5c and 5d portray parallel results at 2 MPa; 
in both simulations the pressure is in the subcritical 
regime. Inspection of the plots shows that for a pure 
mixture at subcritical conditions Le and Le^ff have 
the same values in regions of weak gradients, thereby 
validating our model for Le ^/ / . Quantitative and qual- 
itative differences appear in mixture regions where gra- 
dients are important, and Le// can be as much as a 
factor of 2 larger than Le. With increasing pressure, 
while still in the subcritical regime, these differences 
become larger indicating that the estimate provided by 
Le regarding the relative importance of heat and mass 
diffusional scales deteriorates as the pressure increases. 

Similar plots of Le and Le^f / appear on Fig. 6 (a - 
d) for the low temperature simulations validated above. 
Figures 6a and 6b are for the same pressure as in Figs. 
5c and 5d (2 MPa), but at 445 K, whereas Figs. 6c and 
6d are for 5 MPa (supercritical pressure) and 495 K. For 
the low temperature, subcritical pressure Le and Lccj / 
are in close agreement (Figs. 6a and 6b), again vali- 
dating our model for Lcg / / and comparisons with Figs. 
5c and 5d shows that as expected, the lack of strong 
heating in Figs. 6a and 6b decreases the gradients and 
increases the agreement of Le and Le^f f. However, even 
at low heating, under supercritical pressure conditions 
Le// becomes substantially different from Le as shown 
in Figs. 6c and 6d. Mistakenly, Le indicates that the 
drop is gaseous (Le < 1), whereas Le^fj correctly in- 
dicates that the drop is a dense gas with properties 
approaching those of a Liquid {Le^ff >1). This com- 
parison emphasizes the importance of an appropriate 
model for the determination of the ratio of heat to mass 
diffusional scales which is both of academic and prac- 
tical interest in that design engineers often determine 
combustions chamber dimensions based upon such es- 
timates. 

Fortuitously, the far field boundary conditions of 
the simulations (taken to be those of the data) pre- 
sented in Figs 5 and 6 allowed the uncoupling of the 
pressure and temperature effects on the Le versus 
Le^ff. Wo further explored the combined effect of pres- 
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Conclusions 

A model of fluid behavior under both sub- and super- 
critical thermodynamic conditions has been discussed 
with particular emphasis on the different physics ac- 
cording to the initial conditions with respect to the crit- 
ical point. The model has been exercised for a fluid drop 
for which data are available for model validation. The 
drop is typically colder than its surroundings whose far 
field conditions are prescribed. In the subcritical regime 
and for large emission rates from the drop, there exists 
a film layer in the irmer part of the drop surface and the 
solution of the equations has a convective-diffusive char- 
aicter. In the supercritical regime, there is no material 
surface to foUow, and this introduces an indeterminacy 
in the boundary conditions. To resolve this indetermi- 
nacy one must follow an arbitrary boimdary of interest 
which is here that of the initial fluid drop. The solu- 
tion has then a pure diffusive character, and from this 
solution we calculate the location of the highest density 
greuiient which we identify with the optically observable 
fluid drop. 

The model was exercised for a heptane drop in ni- 
trogen because of the existing data available for com- 
parison. Simulations obtained with this model were 
validated with microgravity experimental data for large 
drops over a wide range of temperatures and pressures. 
The large temperature data were used to determine the 
value of the thermal diffusion factor and further valida- 
tions were conducted with this fixed value. The agree- 
ment between predictions and d^data is excellent at at- 
mospheric pressure and becomes fair at supercritical 
pressure, whereas the rate of regression of the point of 
maximum density gradient is remarkably well predicted 
at all pressures. The numerical predictions show that 
the traditional d^-law is obeyed only in the subcritical 
regime. As the pressure is increased, <P becomes non- 
monotonic with time, with a slope whose magnitude 
increases as a function of time. Thus, we initially iden- 
tify a heating period during which the drop size may 
increase, followed by a period during which the size is 
continuously reduced. The duration of the heat-up pe- 
riod increases with far field pressure. 

Analysis of the results shows that the classical 
Lewis number is not a good indicator of the ratio of heat 
to mass diffusion at high pressure or in the presence of 
strong gradients. Instead, an effectve Lewis number 
previously defined is shown to give a qualitatively con- 
sistent indication of the relative scales of heat and mass 
diffusion. 
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Fig. 1 High temperature comparisons. iZj = 0.35 
mm; /ij = 4 mm, — 0 and = 300 K. In the 
far held Te and p« are specified as in the experiments. 
Simulations at T« = 745 K and O.lMPa, a/ic = 

0.01 ; O.lMPa, olbk = 0.01— o—; O.SMPa, oiik = 

0.01 ; 0.5MPa, as/c = 0.01 - - o - 2MPa . 

Data: 741 K and O.lMPa ■; 749 K and O.SMPa A; 
746K and 2MPa ▼. 



Fig. 2 Intermediary temperature comparisons at 
2MPa. = 0.35 mm; = 4 mm, F® = 0 and 

= 300 K. Simulations at 655 K; ajK = 0.01 : 

0.3 - - - : - 0.3 ; - 0.6 ;0.6 ; agK = 0.01 

— o — . Data at 656 K: ■. 
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0.56 (Mor) 

0..350 

0.90 

0.5 

470 

0.117 (Nom) 


0.135 

1.15 

0.5 
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0.356 (Nom) 
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0.437 (Nom) 



0.89 
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0.138 (Nom) 


0.135 

0.98 
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655 

0.424 (Nom) 





445 

0.097 (Sat) 

0.14 (Sat) 



2.0 

450 

0.096 (Nom) 




2.0 


0.475 (Nom) 


0.360 

0.76 

2.0 

745 

0.4 - 1.7 (Nom) 


0.450 
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5.0 

495 

0.2 - 0.13 (Nom) 





Table 1: Maximum regression rate of the maximum density gradient location, K in mm^/s, obtained from the 
current model (ap), Nomura et al.’s microgravity experimental data (Nom), Sato’s microgravity and normal gravity 
experimental data (Sat), and Morin et al.’s normal-gravity data (Mor), The Nomura et al.’s and Morin et al.’s data 
were provided by the authors, and Sato’s values were read on their graph following the directions given in then- 
paper. In the simulations = 300K and dfi = 0.7mm, while Nomura et al.’s d® was 0.6 - 0.8mm, Sato’s was 1 mm , 
and Morin et al.’s was 1 - 1.5 mm. 




Fig. 3 Intermediary temperature comparisons. 
= 0.35 mm; = 4 mm, = 0 and 7^^ = 300 

K. Simulations at 655 K: O.lMPa ; 0.5MPa ; 

IMPa ; 2MPa . Data; 648 K and O.lMPa ■; 

655 K and 0.5MPa A; 669 K and IMPaT; 656 K and 
2MPa«. 


Fig. 4 Low temperature comparisons. - 0.35 
mm except at 445K where = 0.5 mm; = 4 mm, 
= 0 and = 300 K. Simulations at 470 K; 0.1 

MPa ; 0.5 MPa ; 1 MPa — • — ; at 445 K and 

2 MPa ; at 495 K and 5 MPa . Data: 471 

K and 0.1 MPa ■; 468 K and 0.5 MPa A; 466 K and 1 
MPa ▼; 445 K and 2 MPa 0; 452 K and 2MPa ►; 493 
K and 5 MPa •. 



































































